Conductance and Thermopower of Ballistic Andreev Cavities 
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When coupling a superconductor to a normal conducting region the physical properties of the 
system are highly affected by the superconductor. We will investigate the effect of one or two 
superconductors on the conductance of a ballistic chaotic quantum dot to leading order in the 
total channel number using trajectory based semiclassics. The results show that the effect of one 
superconductor on the conductance is of the order of the number of channels and that the sign of the 
quantum correction from the Drude conductance depends on the particular ratios of the numbers of 
channels of the superconducting and normal conducting leads. In the case of two superconductors 
with the same chemical potential we additionally study how the conductance and the sign of quantum 
corrections are affected by their phase difference. As far as random matrix theory results exist these 
are reproduced by our calculations. Furthermore in the case that the chemical potential of the 
superconductors is the same as that of one of the two normal leads the conductance shows, under 
certain conditions, similar effects as a normal metal-superconductor junction. The semiclassical 
framework is also able to treat the thermopower of chaotic Andreev billiards consisting of one 
chaotic dot, two normal leads and two superconducting islands and shows it to be antisymmetric in 
the phase difference of the superconductors. 



I. INTRODUCTION 

Transport problems have always attracted a lot of 
attention in condensed matter physics. While the 
Landauer-Buttiker formalism which connects the electri- 
cal current to the quantum transmission probabilities of 
a conductor is of key importance for transport through 
nanosystems, similar formulae have also been derived for 
hybrid structures consisting of normal conducting (N) re- 
gions connected to superconductors (S)i^— in which An- 
dreev reflection* plays a crucial role. 

Andreev reflection* can occur whenever a normal 
metal region is coupled to a superconductor. If an elec- 
tron hits the normal metal-superconductor (N-S) inter- 
face with an energy closely above the Fermi energy an 
additional electron-hole pair can be created, and the 
two electrons enter the superconductor forming a Cooper 
pair. The hole however has to compensate the momen- 
tum of the original electron, therefore it retraces the elec- 
tron path. Moreover the hole picks up a phase equal to 
the phase of the macroscopic superconducting wave func- 
tion. 

The early theoretical and experimental investiga- 
tions of transport properties focused on the current 
through the interface of normal metal-superconductor, 
normal metal-insulator-superconductor (N-I-S) and S-N- 
S junctions^. For these the BTK-theorj>i applies, which 
is based on the Landauer type equation 
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where / is the current through the N-S interface with an 
applied voltage V and a measure of the area of the junc- 
tion. In ([1]) i?o is the probability for normal reflection, 
Ra is the probability for Andreev reflection, and / is the 
Fermi function. The BTK theory predicts for N-S junc- 



tions with sufficiently large barrier strengths at the N-S 
interface that the differential conductance d//dV^ van- 
ishes for voltages smaller than the superconducting gap 
A/e. In this regime the conductance is doubled compared 
to the conductance of the same normal conducting re- 
gion with a normal conducting lead instead of the super- 
conducting one: an indication of the proximity eSec^^. 
When increasing the voltage the differential conductance 
has a peak at eV w A and finally approaches the con- 
ductance of the normal conducting region without the su- 
perconductor. However, the total value of the current for 
high voltages exceeds that of a metallic junction by the 
so-called excess current. The early experiments on N-I-S 
junctions were in agreement with BTK-theory. However 
later experiments^»2 found an enhancement of the differ- 
ential conductance at later known as the zero bias 
anomaly. 

Recently Whitney and Jacquod*° considered a some- 
what different type of setup. They considered a bal- 
listic normal conducting region with a boundary giv- 
ing rise to classically chaotic dynamics. Andreev reflec- 
tion and interference between quasiparticles with slightly 
different paths lead to a hard gap in the density of 
states of such chaotic ballistic conductors coupled to a 
superconductor—"—. In Ref. [13 such a chaotic Andreev 
quantum dotis coupled to two normal conducting and 
one superconducting lead and its transport characteris- 
tics was studied. Using a trajectory based semiclassical 
method they calculated the average conductance between 
the two normal leads of such chaotic Andreev billiards up 
to second order in the ratio A^s /An where A's is the total 
number of superconducting channels and N-^ — Ni + N2 
the sum of the number of channels in the normal leads. 
If the superconducting chemical potential is the same as 
that of one of the two normal conducting leads (abbrevi- 
ated to 'superconducting lead' and depicted in Fig. 1(d) ) 
they found that the correction to the classical conduc- 
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FIG. 1. Schema of various Andreev billiards considered here: 
(a) Andreev interferometer with two superconducters, (b) 
double dot setup, (c) chaotic quantum dot coupled to one 
superconducting island, (d) the case of a superconducting 
lead with the same chemical potential as the right lead, (e) 
the so-called "symmetric house", (f) the "asymmetric house" 
where at lead 1 a neck is additionally inserted compared to 
(e). 



tance be huge (of order of iV = TVn + A^s) compared to 
usual weak localization effects, in particular it was shown 
that the quantum correction may become negative or 
positive depending on the ratio Ni/N2. A similar change 
in the sign of the quantum correction to the conductance 
may be caused by a change in the transparencies of the 
leadsi^. 

Using the same approach Whitney and Jacquod 
showed furthermore that to leading order in /A^n the 
thermopower of a chaotically shaped normal metal quan- 
tum dot with two normal leads and two superconduct- 
ing islands (called a "symmetric house" and depicted in 
Fig. 1(e) I is antisymmetric in the phase difference of 
the superconductors. They also argued that the ther- 



mopower vanishes if the two superconductors carry the 
same amount of channels as long as no symmetry break- 
ing neck is inserted at one of the two superconductors 
(c./. Fig. [1(f)). 

Here we combine the tra jectory based semiclassical ap- 
proaches of Refs. 10, 12, 13 and .15 and provide a com- 
prehensive calculation of the conductance and the ther- 
mopower of Andreev billiards. In Refs. [H and [l^ a 
method was developed for the systematic evaluation of 
multiple sums over electron and hole type orbits arising 
in a semiclassical approach to the proximity effect on the 
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FIG. 2. The diagonal diagrams contributing to the conduc- 
tance up to third order in the number A's of channels of the 
superconductor. Here, e and h denote electron type and hole 
type quasiparticles and the asterisk denotes that the path en- 
ters the calculations with the complex conjugated factors. 



density of states of Andreev billiards. Here we further 
extend this recent approach to the conductance. To this 
end a diagonal backbone is introduced which is given by 
a path and its complex conjugated. The quantum correc- 
tion in leading order in 1 /N is then obtained by attaching 
an even number of so called trees (or complex conjugated 
trees) as those used in Ref. [l^ In this diagrammatical 
language, in Ref. [13 the authors restricted themselves to 
at most two trees consisting of just one path pair. There- 
fore their results are valid only for small Ng/Nj^ and the 
validity of their results for larger Ns is not known. Un- 
like the results in Ref. [l6l where the anthers considered 
the distribution of the conductance of chaotic quantum 
dots with one open channel per lead, our results will be 
valid for large numbers of channels in the normal leads. 

We will derive the conductance of the two setups in 
Ref. [13 - na mely the setup with a superconducting is- 
land (see Fig. 1(c) I and with a superconducting lead (see 
Fig. l(d)[ ) - to all orders in Ng/N-^. To this end we start 
llVl bv considering the semiclassical diagrams 



in Sects. Illl- 

and their contribution to the transmission probabilities 
and thus to the conductance to leading order in 1/N. 
In Sect. [V| we apply this approach to the setup with a 
superconducting island. We show that our semiclassical 
result for the conductance coincides with previous ran- 
dom matrix theory results^^ existing for zero magnetic 
field and temperature (though still with a phase differ- 
ence (j) = (f>i — (f>2 between the superconductors). We 
furthermore consider the magnetic field and temperature 



dependence of the conductance of setup Fig. 1(c) 



For the other setup of an Andreev billiard coupled to 



one or two separate superconducting leads (Fig. 1(d) ) we 



will show as a main result in Sect. IVII that the quan- 
tum correction to the classical value of the conductance 
changes its sign not only with the ratio of the number of 
channels in the two normal conducting leads N1/N2 but 
also by tuning the ratio x = Ns/Nj^i. This sign change 
wa not anticipated in Ref. [lol . since it requires an analysis 
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FIG. 3. Pairs of paths contributing to the third-order term in a; = Ns/N-^ of the transmission. Electron (hole) paths are green 
(red). The solid (dashed) lines belong to 7 (7'). A trajectory pair entering from the left and exiting to the right can connect the 
two normal conducting leads while a trajectory pair entering and exiting at the same side can only contribute if the incoming 
and outgoing channel belong both to the same lead. 



to higher orders in x. This conductance correction is also 
shown to oscillate with the phase difference (f> between the 
two superconducting leads with period 2tt. Finally, we 
study the dependence of the conductance on an applied 
magnetic field and temperature. The effects we observe 
for some combinations of the ratios x and N1/N2 turn 
out to be fairly similar to those found in the structures 
containing only one normal conducting lead. 

In Sect. IVIII we show how the methods derived before 
can be extended to calculate the transmission coefficients 
of two dots connected to each other by a neck wehre each 
dot has one further normal and one superconducting lead 
(see Fig. 1(b)). The conductance of this setup is shown 
to also be symmetric in the phase difference. The sign 
of the quantum correction depends on the ratios x and 
n — Nn/ {Ni + N2), where Nn is the channel number of 
the neck. 

In Sect. IVIIII we finally apply our calculations to the 
thermopower of the setup shown in Fig. 1(e) with both 
equal and different numbers of channels as well as to the 



setup shown in Fig. 1(f) We find that for the symmetric 



house with different channel numbers and for the anti- 
symmetric house the thermopower is antisymmetric in 
the phase difference. 



II. CONTRIBUTING DIAGRAMS 

We will evaluate the quantum transmission between 
two normal conducting leads coupled to a classically 
chaotic, ballistic quantum dot which is additionally con- 
nected to superconducting leads such as depicted in 
Fig. [TJ In a trajectory-based semiclassical approach 
the transmission probabilities may be written a.si^— 



b=l a=l C,C' 



where the * denotes complex conjugation. Here, a and 
b label the channels in lead i £ {1,2} and j G {1,2}, 



respectively. ( and C are classical trajectories starting 
at channel a and ending at channel b. The amplitudes 
Aq contain the stability of the trajectory ( and is its 
classical action. Moreover, Th is the Heisenberg time the 
time dual to the mean level spacing. 

We are interested in the conductance averaged over the 
shape of the quantum dot or an energy range small com- 
pared to the Fermi energy but large enough to smooth 
out fluctuations. Moreover we will take the semiclassi- 
cal limit h 0. Therefore the energy dependent action 
difference in ([2]) causes fluctuations cancelled on aver- 
age unless it is of order of h. Thus we have to pair 
the trajectories in such a way that their action differ- 
ence becomes sufficiently small. The easiest way to do 
this is to require that C — C which is known as the di- 
agonal approximation^^. In Fig. [2] the trajectory pairs 
contributing to the diagonal approximation are drawn 
schematically for up to three Andreev reflections. The 
contributions of the diagonal pairs to the conductance 
provide the classical conductance^^ 
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in the case of the superconducting leads. However as 
shown semiclassically in Refs.jl^ for the density of states 
and[l3for the conductance of Andreev quantum dots one 
has to go beyond the diagonal approximation to fully ac- 
count for quantum effects. The non-diagonal trajectory 
pairs contributing to the conductance of normal junctions 
in the limit fi. — have been first considered in Ref. [13 
and generalized to higher orders in 1/N in Ref. Il^: There 
are small regions in which an arbitrary even number - say 
21 - of trajectory stretches come close to each other. I of 
these trajectory stretches 'cross' each other in this region 
while the remaining I ones avoid crossing. The difference 
between a trajectory ^ and its partner C then leads to 
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FIG. 4. Diagrams corresponding to the trajectory pairs shown in Fig. |3] The full circles denote encounters while the empty 
circles denote Andreev reflections. Note that an encounter touching the superconductor is marked as Andreev reflection. An 
encounter touching a normal conducting lead is shown as an empty box. The solid (dashed) line represents {('). 



a small action difference as long as these stretches are 
close enough to each other. Such a region with I cross- 
ing trajectory stretches and I trajectory stretches 'avoid- 
ing crossing' will be referred to as ^-encounter. Between 
these ^-encounters two different trajectory stretches re- 
trace each other forming a path pair with vanishing ac- 
tion difference which will also be referred to as a link. 

In what follows we will identify the relevant trajectory 
pairs contributing to the conductance beyond the diago- 
nal approximation in leading order in the inverse channel 
number 1/N. The diagrams with two Andreev reflections 
may be found in Ref. [lO|- However we want to go beyond 
second order in a; = Ng/N^- The trajectories contribut- 
ing in third order in x, i.e. trajectories with three An- 
dreev reflections, are shown in Fig. |31 The first task is to 
find a structure in the diagrams contributing at leading 
order in the channel number. To facilitate this we can 
redraw our semiclassical diagrams in a skeleton form and 
represent encounters and path pairs by nodes and lines. 
For example the diagrams contributing to third order in 
X, shown in Fig. |31 can be redrawn as in Fig. 21 

We first consider how to read of the channel number 
dependence from a given diagram, i.e. we use the dia- 
grammatic rules used in Ref. '23' disregarding an energy 
and magnetic field dependence and any signs for the mo- 
ment. A path pair hitting lead j contributes a factor of 
channel number Nj. The path pair, or link, itself how- 
ever contributes a factor 1 /N while each encounter con- 
tributes a factor N. From the trajectory pairs shown in 
Fig. [3] we see that if we cut off all e-h and e*-h* pairs 
we again get a diagonal like contribution as depicted in 
Fig. [5l For example if we cut the e-h pair at the very left 
of ee3I we get the diagonal contribution to the second or- 
der in X in Fig. [5] since there are two Andreev reflections 
and if we cut the 'off-diagonal' parts in say eeSIII we get 
a diagonal contribution at first order in x. 

Staying at leading order in 1/N implies that the 'off- 
diagonal' path pairs cannot consist of one ^- and one 
stretch since each of those path pairs has to be traversed 
by ^ and C in the same direction. Thus in order to come 
back from the off-diagonal part starting with an C-C pair 
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FIG. 5. If the e-h path pairs are cut off a diagonal type 
diagram remains. 




FIG. 6. Diagrams we neglect in leading order in due to 
the formation of loops: (a) A non-diagonal (-C^' pair causes 
the formation of a loop, (b) A loop formed by an off-diagonal 
e-h path pair, (c) A loop formed due to the lack of a diagonal- 
type (-(^' path pair. 



we have to connect this off-diagonal part to the diagonal 
'backbone' by a sec ond ^- C' path pair thus forming a loop 

This loop however adds a link. 



as indicated in Fig. 6 (a 



thus giving a factor 1/N compared to the same diagram 
without the loop, and therefore decreases the number 
of Andreev reflections by at least one and therefore the 
contribution to the conductance is suppressed by a factor 
of the order 1 /Ng such that it would contribute to sub- 
leading order in the inverse channel number. Therefore 
the 'off-diagonal' parts may only consist of e-h or e*- 
h* pairs. In the same way we may neglect loops formed 
by e-h or e*-h* path pairs as the one in Fig. 6(b) 



In terms of graphs, the 'off-diagonal' parts again be- 
come rooted plane trees as in Refs. [H M- These trees 
start with a link (root) which connecting an encounter to 
the diagonal like backbone. From this encounter several 
further links emerge all ending again at an encounter or 
at a superconducting channel. In contrast to the trees 
in Refs. [l^ IHj the trees here - we will call them 'side- 
trees' - start at the 'diagonal encounter' such that their 
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root does not touch a channel but the diagonal backbone 
instead. Note that we draw the diagrams such that the 
non-complex conjugated side trees are at the upper side 
of the diagonal backbone while the complex conjugated 
ones are on the lower side of the diagonal backbone. 

The fact that the path pairs along the backbone are 
composed only of (^-(' pairs is again due to neglecting 
loops: The two trajectories C and C must both start at 
lead j and end at lead i. Therefore the path pairs hitting 
the normal leads have to be pairs and so, if there is a 
'diagonal' encounter entered by a C,-C' pair and left only 
by e-h and e*-h* pairs, there must be a corresponding 
encounter entered only by e-h and e*-h* pairs and left by 
a pair. Therefore we again would get a loop essen- 
tially formed by one e-h and one e*-h* pair as shown in 



Fig. 6(c 



All told the diagrams have to consist of a diagonal type 
'backbone' consisting of C-C path pairs and encounters 
(which may also touch the superconductor) and C- and C'- 
side trees emerging from these diagonal encounters. Note 
that when pairing a C with a C stretch these stretches 
have to be traversed by the same kind of quasiparticle, 
i.e. it has to be an e-e* or a h-h* pair. This is related 
to the fact that each encounter has an even number of 
entering and exiting path pairs. 

However, there is still one possibility left we have not 
mentioned yet but that needs a special treatment. If 
the 'diagonal' part consists of only two path pairs and 
one 2-encounter with one C-side tree (a side tree formed 
by C) and one C'-side tree this encounter can be moved 
into one of the normal conducting leads, say lead i. An 
example for an 2-encounter touching the incoming lead is 
the trajectory labelled by heSIV in figures [3] and S] which 
arises from the trajectory labelled by heSVII by moving 
the encounter into the lead. However this is only possible 
if the trajectory connects lead i to itself and thus if the 
electron is coherently back scattered. In this case we have 
only one side tree and one complex conjugated side tree 
but no 'diagonal' part. 

Since we know the structure of the trajectory pairs con- 
tributing at leading order in the inverse channel number 
we can start evaluating them. Because the contributions 
of the encounters and the stretches are multiplicative^- 
we may factorise the contribution of a given diagram into 
the contributions of side trees starting at the first en- 
counter with an a- type quasiparticle, P°'{e,x), the first 
encounter itself and the diagram remaining when cutting 
the diagram after the first encounter as in Fig. [51 We 
will first evaluate the contribution arising from the sum- 
mation over a possible side tree. 



III. SIDE TREE CONTRIBUTIONS 

We restrict ourself to sufficiently low temperatures 
such that only energies cEt (measured with respect to 
the Fermi energy and in units of the Thouless energy 
Et = h/2TD where td is the mean dwell time) much 



smaller than the superconducting gap A have to be taken 
into account eE^ <g: A. This allows us to approximate 
exp[— iarccos(e/A)] « — i such that the scattering ma- 
trix of Andreev reflection becomes independent of the 
energySi. Thus the diagrammatic rules for the C-side 
trees readi^^^^ 

• An e-h path pair contributes [A^ (l + ie + fe^)] ^ 

• An /-encounter contributes — A^ (l + ile + Pfe^) 

• An e-h path pair hitting the superconductor Sj con- 
tributes A"q . . 

• An /-encounter touching the superconductor Sj 
contributes A"Sj . 

• Each Andreev reflection at the superconducting 
lead j converting an electron into a hole contributes 

• Each Andreev reflection at the superconducting 
lead j converting a hole into an electron contributes 
-ie'"*^ . 

with b (X B/h where B is the magnetic field applied 
to the system. The proportionality factor depends on 
the actual system^'^. These diagrammatic rules have to 
be complex conjugated for a C'-side tree and imply that 
when exchanging electrons and holes we just have to re- 
place 4> ^ —4>- Thus a side tree starting with a hole 
gives the same contribution of a side tree starting with 
an electron but with negative phase. Therefore we only 
need to evaluate side trees starting with electrons. 

The evaluation of the side trees then follows essen- 
tially those in Refs. [l^ [l^ and [2^ However here the 
root of the tree does not hit any channel and therefore 
can not touch the superconductor which simplifies the 
calculation. Moreover from the rules above for a path 
pair hitting a channel in the superconductor Si we get 
a factor — ie~"^/^A"si if an electron hits the channel and 
_igi0/2 jYg^ if a hole hits the channel, rather than just a 
factor of the numbers of channels, and equivalently for a 
path pair hitting S2. 

Similar to Ref. [ij as long as the phase difference 
is zero and no encounter touches the superconductor the 
contribution of a side tree with characteristic v - which is 
the vector which /-th entry is the number of /-encounters 
of the tree - but without the contribution of the path 
pairs hitting one of the superconductors is 



„=i (1 + ie + fo^ 



2i„-l 



V'(v) 



(i+i.+6^)-"n 



a=i ll + ie + b' 



(5) 



Here the encounters have been labelled by a and we used 
that the side tree has to satisfjsi^ n = i(v) — V{v) + 1, 



where n is the number of hnks touching the supercon- 
ductor, V{v) = ^ is the total number of encounters 

l>2 

of the tree and L(v) — 'Y^lvi. This is because every l- 

l>2 

encounter creates 21 — 1 additional path pairs and each 
path pair has to end either in an encounter or at the 
superconductor. 

We then enumerate the number of Z-encounters by xi 
and the number of Z-encounters touching the supercon- 
ductor Si at an odd numbered channel by z^'j. An l- 
encountcr touching the superconductor means that the I 
incoming trajectory stretches hit the superconductor at 
one and the same channel. We look at the generating 
function /(x, Zo"'^'' , zi^'') which counts the number of pos- 
sible side trees and their encounter types and derive a 
recursion relation for it by cutting the side tree at its 
top encounter. If the top encounter is traversed by 21 
stretches and does not touch the superconductor the tree 
then has the contribution of the top encounter times that 
of all 21 — 1 subtrees giving xif^ , where / is the same 
as / but with replaced by —0 accounting for the fact 
that each even numbered subtree starts with a hole in- 
stead of an electron. If the top encounter however is an 
encounter traversed by 21 stretches touching its con- 
tribution is Zq*;/'^^. In total we therefore have 
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(6a) 



(6b) 



where the first two terms account for empty side trees 
which consist of one link and one Andreev reflection at 



Si or S2 and z'^\ is the same as but with (j) replaced 
by -(j). 

Due to the fact that the links of the side trees are tra- 
versed by one electron at energy -\-eh/2TD and one hole 
at energy — eft/2rD in opposite directions an /-encounter 
consists of I electron-stretches traversing the encounter 
in the same direction and I hole-stretches traversing 
the encounter in the opposite direction. Thus we have 
xi = - (1 + i/e + Z^fo^) / (1 + ie + 5^)' f'-i in line with 
([5]) . The powers of f are included in order to keep track 
of the order of the trees. Now consider an /-encounter 
touching Si. According to the diagrammatic rules after 
extracting the factor (1 -I- ie -I- 6^)~" as in (O the con- 
tribution of the encounter and the link connecting the 
top encounter to the backbone is Ng-^/N. However we 
have to include the phase factors contributed by the An- 
dreev reflections. To evaluate this phase factor we look 
at the /-encounter touching the superconductor as aris- 
ing from an /-encounter inside the dot by sliding it into 




FIG. 7. If the odd numbered subtrees have zero characteristic 
and hit the same superconductor the top encounter may be 
slid into the superconductor. 



the superconductor as indicated in Fig. [71 This is only 
possible if the odd numbered subtrees have zero charac- 
teristic and hit the same superconductor. By sliding the 
encounter into the superconductor the total number of 
Andreev reflections do not change such that the phase 
factor provided by the encounter touching the supercon- 
ductor is the same as the phase factor provided by the 
odd numbered side trees we start from. For a side tree 
starting with an electron the odd numbered subtrees with 
zero characteristic provide one Andreev reflection con- 
verting an electron into a hole. Since from an /-encounter 
/ odd numbered subtrees emerge the phase factor of an 
/-encounter touching Si is — ie"''"^'. 



^(_i)^_^e 



'ii0/2~;-l 



J2) 



(_i)'^e"^/2~;-i 



N 



(7a) 
(7b) 



The total power of a tree with 2n — 1 Andreev reflec- 
tions is again — l)vi = L — V = n — 1. Thus in order 
to get the required prefactor of (1 + ie + fo'^)"" we can 
make the change of variables 



/ = g(l + ie + 62), 



1 



(8) 



After making this change of variables and performing 
the summations in ([Sai b) using geometric series we get 
in view of Eqs. (ITal b) 

(l + i£-Kb")g ^ {2b^ + ie)rgg^ ^ {1 + rgg) rgg^ 
1 - rgg (l „ rggf 



+ 



ix{l + y) 



ix(l 



- rgg) 
V) 



2(1 + x) (e'"^ 



irg 



2(1 + x) (e 



i0/2 , 



\rg 



(9) 



and the same equation with g and g exchanged and re- 
placed by —cf). Here we used N^/N^ — x and introduced 
the difference of the numbers of channels of the two su- 
perconductors y = (iVsi — Ng^)/Ns such that y = cor- 
responds to the case of equal numbers of channels and 
y = ±1 to the case of just one superconductor. 
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In the case that the two superconductors provide the 
same number of channels (y=0) those two equations are 
the same implying g — g and ([9]) is equivalent to an 
algebraic equation of 7th order in g. This increase in the 
order of the equation with respect to the same case for 
the density of states^^ is due to the fact that in the case 
of the density of states we had no normal leads. 

The contribution P'^ of the side trees starting with an 
electron is then obtained by giving all trees the same 
weight by setting r — I'm g. The contribution of the side 
trees starting with a hole are then given by replacing (p 
by —0 in (; or setting r = 1 in ^. After setting r = 1 
and eliminating say P'' the contribution of the side trees 
starting with an electron P"^ is given by a rather lengthy 
equation of in general 11th order which factorises in the 
case y = such that P*^ = — P \s given by 

-P^ + (2i/3 + ifix) P^ 

+ {-b^x + 3 + iea; - 5^ + ie) P^ 

+ {-i(3x + 216^/3 + 2e/3 + 2ih^l3x + 2e/3a; - Aip) P* 

+ (-2ie - 3 - 2iea;) P^ 

+ (216^/3 + 2i52/3x - 2el3x - 2e/3 + 2\I3 - i/3a;) P^ 

+ {lex + 6^ + + 1 + ie) P + i(ix = 0. 

(10) 

If no magnetic field is applied (5 = 0) the equation may 
be factorised, and one has to solve an equation whose 
order is lowered by 2. 

If the Andreev interferometer consists of two supercon- 
ductors with the same numbers of channels (y=0) the side 
tree contributions only depend on /? = cos((/)/2) rather 
than on (j) itself. Therefore in this case the side tree con- 
tributions are symmetric in (j) and the contribution of a 
side tree starting with a hole is the same as that of a side 
tree starting with an electron. In the most simple case 
of the absence of a magnetic field, zero temperature {i.e. 
e = 0) and zero phase difference reduces to a second 
order equation: 

- P"^ + iP + \Px - x (11) 

yielding 

P(0, x)^]^[l + x-y/l + Qx + x'^Y (12) 

Note that we take the solution satisfying P(0,0) — 
since when there is no superconductor the correction of 
leading order in the channel number has to be zero. 

IV. TRANSMISSION COEFFICIENTS 

We will now demonstrate how to calculate the trans- 
mission coefficients T^"'' for transmission from lead j to 
lead i while converting an a-type quasiparticle in a {3- 
type one, using Tf^ as an example, as the evaluation of 
the other transmission coefficients will be similar. 




FIG. 8. (a) A diagram contributing to Xi is split right af- 
ter the first 4-encounter and decomposes into two separate 
diagrams where the second one contributes to Ti-^. (b) To 
sum over all diagrams starting with an Z-encounter we can 
remove a factor corresponding to the first encounter (and its 
side trees) and a sum again over the transmission diagrams. 



We first order the sum over all diagrams contribut- 
ing in leading order in the channel number with respect 
to the first encounter. Then the first summand is of 
course the diagram corresponding to the upper left tra- 
jectory in Fig. [5] Next there are all diagrams whose 
first encounter is a 2-encounter followed by all diagrams 
whose first encounter is a 3-encounter etc. Note that 
we also allow for the first encounter to touch the su- 
perconductor or (if the first encounter is a 2-encounter 
and i = j) the normal lead. We denote the contribution 
of the sum over all diagrams having an ^-encounter as 
their first encounter and contributing to T^^^ in leading 
order in the number of channels by A/. We may include 
the diagonal diagram without any encounter by setting 
Xi ~ NiNj/N. The transmission coefficients are then 
given by T^J' — J2i>i-^i- Now we fix Z > 2 and split 
all diagrams contributing to A; right after the first en- 
counter into one part consisting of the first path pair and 
the first encounter together with its side trees and the re- 
maining part such as indicated in Fig. |8(b)[ Note that 
the diagonal type path pair leaving the first encounter 
is completely included in the second part. Since the di- 
agrammatic rules are multiplicative the contribution of 
a diagram is given by the product of the two parts and 
hence they all have a common factor which is given by 
the first diagonal type link, the first encounter and the 
side trees emerging from it. To sum over all diagrams 
starting with an Z-encounter we pull out this factor and 
are left with a sum over the transmission diagrams as 
depicted in Fig. |8(b)[ This sum runs over all possible 
diagrams contributing to Tf^^ if the first encounter is left 
by an electron and to TfJ^ if it is left by a hole. How- 
ever in order to be able to fully identify the sum over the 
second parts as the transmission we have to reassign the 
contributed number of channels A^^ contributed by the 
first path pair leaving lead j to the second part. Due to 
the two possibilities of which transmission coefficient the 
remaining diagrams contribute to we split A; into two 
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FIG. 9. Simple examples for encounters touching a supercon- 
ductor. The electron paths are shown green while the hole 
paths are shown red. The solid lines belong to while the 
dashed ones belong to If the quasiparticles entering an en- 
counter touching the superconductor following an diagonal- 
type path pair the diagonal-type path pairs leaving it are 
traversed by holes and vice versa. 

parts 

Xi = AfTif + BfTf^, 

where Ai is the contribution of the first e-e* pair and 
the /-encounter the path pair enters together with all 
side trees and with the entering and exiting quasiparticle 
being the same. Bf is the same but with the entering 
and exiting quasiparticle being different. 

The transmission coefficients may therefore be written 

as 

N-N- °° °° 

T!j-^+Y. + E Bin;", (isa) 

oo oo 

= ^ 7;f -t- J2 Bj'T^- (13b) 

1=2 1=2 

A'l and Bj^ are the same as and Bf , respectively, but 
with electrons and holes exchanged. Equation (I13bp is 
obtained in the same way as (|13ap but with the additional 
condition that there is no diagram without any Andreev 
reflection contributing to it since converting an electron 
to a hole requires at least one Andreev reflection and 
therefore one encounter. The formulae for Tj]^ and T']'^ 
are the same but with e and h exchanged. 

The next task is to find out what causes the encounter 
which is entered by an electron to be left by an electron or 
a hole. The trajectories in Fig. [3] and their corresponding 
diagrams in Fig. |4] indicate that, as long as the first en- 
counter does not touch the superconductor, an encounter 
entered by an electron is left by an electron if the num- 
ber of side trees on each side of the diagonal backbone 
emerging from this encounter is even (such as in the dia- 
grams eeSII and heSV) and by a hole if it is odd (such as 
in the diagrams eeSIII and he3 VII) . If the first encounter 
however touches the superconductor the encounter is al- 
ways left by a hole if it was entered by an electron. This 
is also indicated in Fig. |9] 



We will now show that this indeed holds for all 'di- 
agonal encounters' entered by a diagonal type e-e* pair 
by starting with considering encounters not touching the 
superconductor. Since an Z-encounter connects 21 links 
to each other, each diagonal /-encounter, where 2 of the 
links belong to the backbone, provides in total [21 — 2) 
side trees implying that if the number of ^-side trees is 
even the number of C'-side trees is even too, or they are 
both odd. Furthermore each side tree provides an odd 
number of Andreev reflection and therefore a conversion 
of an electron into a hole or vice versa, since each of its 
/-encounters is left by (2/ — 1) additional path pairs and 
each path pair increases the number of Andreev reflec- 
tions by one (this is closely related to the fact that we 
consider diagrams contributing at leading order in the 
number of channels). Thus, as long as the first encounter 
does not touch the superconductor, the entering electron 
leaves the encounter as an electron if the number of side 
trees p built by C, is even and as a hole if the number of 
side trees built by C is odd. 

However if the first diagonal /-encounter touches the 
superconductor the first side tree starts with a hole in- 
stead of an electron and is therefore left by an electron. 
Since the electron leaving the first side tree hits the su- 
perconductor the second side tree again starts with a 
hole. If one proceeds inductively one finds that every 
side tree starts with a hole and is left by an electron 
which after that undergoes again an Andreev reflection. 
Therefore if the first encounter entered by an electron 
touches the superconductor it is always left by a hole 
and we can view it as arising from an /-encounter with 
an odd number p of C-side trees slid into the supercon- 
ductor as indicated in Fig. [TU] and therefore contributes 
to Bf . An /-encounter may then touch the superconduc- 
tor if the number of ^-side trees p is odd and the odd 
numbered (^-side trees, which are the side trees traversed 
by C, after an odd number of traversals of the encounter, 
as well as the odd numbered C'-side trees have zero char- 
acteristic {i.e. consist of just one link and one Andreev 
reflection) . Moreover the links of the odd numbered side 
trees have to hit the same superconductor such that the 
channels can coincide. When sliding such an encounter 
into the superconductor the channels at which the odd 
numbered side trees hit the superconductor coincide and 
the links vanish. Therefore beside the diagonal-type path 
pairs from such a diagonal /-encounter touching the su- 
perconductor p = {p — l)/2 even numbered ^-side trees 
starting with a hole and [(2/ — 2— p) — l\/2~l~2~p 
even numbered C'-side trees, which also start with a hole 
emerge. 

Thus if we denote the contribution of the first a-a* 
pair and of the first /-encounter inside the dot with p 
C-side trees by x"- and the contribution of the Andreev 
reflections provided by the first /-encounter touching the 
superconductor Sj created by sliding an /-encounter with 
originally p C-side trees into the superconductor Sj by 
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FIG. 10. (a) An 3-encounter may touch the superconductor Si if the odd numbered side trees have zero characteristic and hit 
the same superconductor. The number of Andreev reflections stays the same. If the encounter touches the superconductor an 
entering electron is converted into a hole, (c) A more complicated diagram with two diagonal encounters that may touch the 
superconductor. Note that additionally the fourth side tree may also touch the superconductor but this does not affect the 
diagonal encounter but is instead included in the side tree recursion. 



l-l 



p=0 
1-2 



pn 



p=0 



1-2- 



pn 



p=0 
1-2 



ph 



l-l-p 



p=0 



[{P' 



{pr^' (pr \{p') 1 {inT'-' + e -'2^+1, (pr [ipr] 



,e\*l l-i-p 



l-2-p 



(14a) 
(14b) 
(14c) 
(14d) 



where we have used that p has to be even for Af and 
thus replaced p = 2p and odd for Bf with p = 2p+ 1. 

The next and final step is to find the contribution of 
the encounters. For that we would like to recall the dia- 
grammatic rule for an /-encounter traversed by trajecto- 
ries with e nerg ies ±e and in presence of a magnetic field 
b from Ref. US 

• An Z-encounter inside the dot contributes a factor 
-N{l + -nie + fi^b^). 

Here 77 is the difference between the number of traver- 
sals of e-stretches and the number of traversals of e*- 
stretches and fi is the difference between the number of 
C-stretches traversed in a certain direction and the num- 
ber of C'-stretches traversed in the same direction. Since 
every electron path of the side tree is retraced by a hole 
every second stretch connected to a C-side tree is an e- 
stretch and they are all traversed in the same direction we 
choose arbitrarily as 'positive'. Therefore if the number 
of C-side trees is even the number of e-stretches traversed 
in positive direction is simply p/2. If p is odd we have to 
account for the fact that the first ^-side tree starts with 



an electron and the last one also does. Thus there are 
{p+l)/2 e-stretches traversed in positive direction in the 
encounter. In the same way one finds that the number 
of e*-stretches are {21 -2-p)/2 and [(2/ ~2-p) + l]/2, 
respectively. Since the diagonal path pair is traversed 
by C and C' in the same direction the directions of the 
e*-paths are also positive. Since the holes retrace the 
electron paths their directions in the encounters is nega- 
tive. Thus in both cases one finds that rj = fi = [p—l + l). 
So 



xi^p = - 1 + {p - I + l)ie + {p - I + if 



(15) 



For the contribution z^^ ■ which arises by sliding an l- 
encounter into the superconductor—, as shown in Fig.fTOl 
we remember that the number p of ^-side trees emerg- 
ing from it is odd and the odd numbered C-side trees as 
well as the odd numbered C'-side trees consist of only 
one path pair and one Andreev reflection (i.e. they have 
zero characteristic). Moreover the Andreev reflections of 
the odd numbered side trees have to be all at the same 
superconductor. The contribution of the encounter it- 
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self and the first path pair is then Ns^/N. However we 
also include the factors contributed by the Andreev re- 
flections in zfp J, too, which are stated in Sect. IIIII As 
for the side trees, these phase factors may be determined 
by looking at the odd numbered side trees before sliding 
the encounter into the superconductor since the number 
of Andreev reflections of the C,- and C'-trajectory can not 
change when sliding the encounter into the superconduc- 
tor. Consider the p + 1 — {p + l)/2 odd numbered side 
trees which have zero characteristic and hit say S^: The 
Andreev reflections provided by these side trees convert 
an electron into a hole and thus the Andreev reflections 



each provide a factor — ie~"^' . Hence in total the Andreev 
reflections of the odd numbered C-side trees provide a 
factor (— 1)''+^ e^''^^^-''^' . Analogously the Andreev re- 
flections of the odd numbered C'-side trees contribute a 
factor {'■-P-^e^(.''-P-'^)^i , Thus in the case of two super- 
conductors with phases (pi = —02 = 0/2 the phase factor 

„-i(2p-Z+2)0/2 



included in zfp ^ is given by (— i) 
We thus have 



p[i- 



N^.i- 
N ^ 



-ig-i(p-;+i)0/2^ 



(16) 



For zf p 2 '^6 have to exchange 



by ■ Moreover z?*, 



^i.,p,j\4>->--'t> 



^ and replace A^Si 
-A,. Therefore we have 



i-i 



i-p-i 



- ^ [l + i (2p - ; + 1) e + (2p - ; + 1)' b'] {P- f {P^Y [{Py] '-"^^ [{P'^) 

1-2 

E 

P=0 L 

Xil+ y)e-'(2p-'+2)0/2 ^_iphY ^phy^ '-^-^ ,^(1 _ ^)^i(2p-;+2)0/2 ^_-^phY (p/.) '^"^^ 



(l + \{2p-l + 2)e+{2p-l + 2f b^) {P^f+' {P^Y [{P^Y] [{P'^) * 



(17a) 



2{l + x) 



2(l + a;) 



(17b) 



where we again used y = (A^Si — A^S2)/-^s- The case 
Ns^ = Ng^ is then obtained by setting y — while 
the case of just one superconducting lead corresponds 
to y — ±1. Since exchanging electrons and holes corre- 
sponds to replacing (j> by —0, A'^ and i3f are obtained 
by the same formulae but with (f> replaced by —0 includ- 
ing an exchange P"^ -f-)- P''. The sums may be performed 
using geometric series and yield our main result. Along 
with (|13a[) and (jl3bp it contains all rhe diagrams, and 
their semiclassical contributions, generated recursively. 
Note that if the numbers of channels of the supercon- 
ducting leads are equal the symmetry of P towards the 
phase implies that A" and B" are symmetric in yield 
ing A 



- and Bf 

rpeh 



B^ and thus T'; 



T'f and 



rphe 
ij 

Therefore we now have all the necessary utilities to 
calculate the conductance of Andreev billiards with two 
superconducting islands. When the incoming and out- 
going lead are the same, i = j, and the first encounter 
is an 2-encounter this encounter may enter the lead. In 
this case however the encounter simply contributes Ni 
and the diagrams consist of one C-side tree and one C'- 
side tree. The contribution of these diagrams is therefore 
simply 



The transmission coefficients necessary for calculating 
the conductance may be calculated by evaluating the 
side tree contribution by solving ([9]), inserting this into 
P7a[ b) and performing the sums and finally inserting 
into (llBai b) and solving for the transmission coefficients. 



V. CONDUCTANCE WITH 
SUPERCONDUCTING ISLANDS 



We now evaluate the conductance of Andreev billiards 
with two normal leads. We first consider a chaotic quan- 
tum dot coupled to two normal conducting leads and 
one or two isolated superconductors with equal number 

The chemical poten- 



of channels as shown in Fig. 1(c) 



tial of the superconducting lead is then adjusted by the 
dot such that the net current in the superconductor van- 
ishes. The dimensionless conductance g — irhl/le^V) 
with / the current and V the voltage drop between the 
two normal leads, in this case is given at zero tempera- 
ture hy^ 



SijNi\P^\'^ if the dot is entered by an electron, (18a) 
(5„7Vi|P''p if the dot is entered by a hole. (18b) 



21 



rphe 
^21 



rpherphe rpherphe 
11 ^22 ^21 12 



^ rphe 
^11 



Tihe I rphe 
21 ~^ ^12 



(19) 
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FIG. 11. Conductance correction scaled by A''i as a function 
of the number of channels of the isolated superconductor x = 
iVs/iVN for N2/N1 = 0.2 (dotted line), N2/N1 = 1 (solid line) 
and iV2/iVi = 100 (dashed line) 



A. Low temperature 



(a) 



(b) 



1. One superconductor 



Using Eqs. (|18ai b) in the simplest case without phase 
difference the random matrix result for the conductance 
correction 6 



g — Qcir^ can be reproduced: 



5g^ 



N1N2 



+ i (1 + x)^ - i (1 + x) Vl + 6a; + a;2 



(20) 

The conductance correction is shown in Fig.[TT]as a func- 
tion oi X = Ns /TVn . From ((20)) it can be easily seen that 
the conductance in this case is symmetric in exchanging 
A'^i and 7V2, as one would expect due to the symmetry of 
the setup. Moreover with this setup the superconductor 
always increases the conductance. In the limit of large 
numbers of superconducting channels we find that the 
conductance is doubled compared to the classical limit 
X = and hence approaches the conductance of an N-S 
interface^. 



2. Two superconductors with phase difference 

In Ref. using RMT for a finite phase difference the 
authors could calculate the transmission only numerically 
but for all N. Moreover they restricted themselves to 
the case A^Si = • With our semiclassical approach 
however we are able to calculate it at zero temperature 
analytically for all cases, as long as N is large. Using this 
we could reproduce the large- limit in Ref. [l3 shown 
in Fig. [12] for the case iVg^ — Ns2- Moreover Fig. [T3l 
shows the dependence of the conductance on the channel 
number difference of the superconductors. The conduc- 
tance correction vanishes for tt since the phase ac- 
cumulated at each Andreev reflection causes destructive 
interference. In this case the conductance is a monotonic 
function of the phase difference ^ up to tt. 

The symmetry found in Figs. [12] and [13] results from 
the fact that electrons and holes contribute symmetri- 
cally to the conductance. The 27r-periodicity may also 



FIG. 12. The conductance correction as a function of the 
phase difference for iVsi = A^Sa (a) N2/N1 = 0.2 and x = 0.5 
(space dashed line), x = 1 (dotted line), x = 2 (solid line), 
a; = 3 (dashed line), x — 5 (dashed dotted line), (b) x — 1 
and N2/N1 = 0.2 (space dashed line), N2/N1 = 0.5 (dotted 
line), iV2/A^i = 1 (solid line), N2/N1 = 3 (dashed line) and 
A^2/iVi = 7 (dashed dotted line). 
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FIG. 13. Dependence on the difference of the numbers of 
superconducting channels with y = 1 (space dashed line), 
y = 0.95 (dotted line), y — 0.9 (solid line), y = 0.5 (dashed 
line) and y = (dashed dotted line) for (a) iV2/iVi = 0.2, 
X = 0.5 and (b) N2/N1 = 7, x = 0.2. 



be found using (jl7al b): If we increase (jj by 2tt the side 
tree contribution changes its sign. This does not affect 
Af or the first part of . If Hs odd in the last two terms 
in changing the sign of the side tree contribution re- 
sults in a change of the sign of these two parts. However 
increasing the phase by 2tt also yields an exchange of the 
sign of the phase factors cancelling the change of sign of 
the side tree contributions. If I is even we again have 
an even number of side tree contributions and the phase 
factors also do not change their sign. 

The crossover from two superconductors to just one 
superconductor is smooth and monotonic as shown in 
Fig. [13] We found that the bigger the difference in the 
numbers of channels the faster the amplitude changes. 



3. Magnetic field 

Whitney and Jacquod found in Ref. |l3 that for small 
X the conductance of an Andreev quantum dot with an 
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FIG. 14. The conductance correction for the setup with one 
superconducting island as a function of the magnetic field 
dependence for N2/N1 = 0.2 and x = 0.5 (space dotted line), 
X = 1 (dotted line), x = 2 (solid line), x = 3 (dashed line) 
and X = 5 (dashed dotted line) 



isolated superconductor decays at T = with increasing 
magnetic field as (1 + 6^)"^. For higher orders in x the 
(1 + 6'^)~^-decay mixes up with terms decaying as (1 + 
6^)" with n > 2. This leads to the behaviour shown in 
Fig. [131 It can be seen that the conductance correction 
decays very quickly. Since the magnetic field enters the 
transmission coefficients and therefore the conductance 
quadratically the conductance is symmetric in reversing 
the magnetic field. 



B. Temperature dependence 

If we want to allow for non-zero temperature each 
transmission and reflection coefficient in (|19p has to be 
multiplied by the negative derivative of the Fermi func- 
tion and integrated over energy. We evaluated these inte- 
grals numerically using Gaussian quadrature with a to- 
tal accuracy 10~^° and a truncation of the integral at 
e — 1009 with 9 — kBT/ET being the temperature mea- 
sured in units of the temperature corresponding to the 
Thouless energy. Doing so we flnd that the superconduct- 
ing island obeys a monotonic temperature dependence: 
The conductance correction has its maximum at T = 0. 
For higher temperature it is damped due to the mixing 
with higher energies for which the side tree contributions 
become smaller because of the loss of coherence of the 
electrons and holes. As the temperature tends to infinity 
the conductance correction vanishes slowly. In Fig. [T5] 
we plotted the conductance correction of the setup with 
one superconducting island versus the temperature. 



VI. CONDUCTANCE WITH 
SUPERCONDUCTING LEADS 

Next we consider superconductors with externally con- 
trolled chemical potential. In particular we will consider 
superconductors lying on the same chemical potential as 
one of the two normal conducting leads, say lead 2. Such 
a setup is schematically shown for the case of one super- 



FIG. 15. Temperature dependence of the conductance cor- 
rection for the setup with an superconducting island for 
iV2/A^i = 0.2 and a; = 0.1 (dotted line), x = 1 (solid line) 
and X = 2 (dashed line). 



conducting lead in Fig. 1(d) The current in lead i may 
be calculated by the Landauer-type expression2i2£ 



df 



de 



(21) 



where / = [exp(— e/0) -t- 1] is the Fermi function with 
the temperature again measured in units of the Thou- 
less energy E^, 11 j and 113 are the chemical potentials in 
the normal conducting lead j and in the superconductor, 
respectively. 

Of course we could use the transmission coefficients 
themselves calculated in Sect. IIVI However we would 
like to present a slightly different way to calculate the 
conductance here which in the case of the superconduct- 
ing leads simplifies the calculation. For simplicity we 
will present this way only for the case that the numbers 
of channels of the superconducting leads are equal since 
the modifications one would have to do in order to in- 
clude different numbers of channels are the same as in 
Sect. IIVI According to (|21l) we have to calculate the dif- 
ference between the Andreev and normal transmission, 

namely ff^ = T'^f - Tlf and f^) = Tf^ - T'f. To do this 

*j *j ''J ''J ''J 

(considering only the case of an incident electron) we es- 
sentially perform the same steps as before and split the 
diagrams at their first ^-encounter. In the same way as 
above we find that the sum over the remaining diagrams 
again contributes to Tf^ or 3^^ depending on whether the 
quasiparticle leaving the encounter is an electron or a 
hole. An additional sign arises as follows: Consider for 
example a diagram contributing to T^-" thus contribut- 
ing to with positive sign. Then if the number of 
(■-side trees arising from the first encounter is even the 
remaining diagrams contribute to T^^, too, and therefore 

it again contributes with a positive sign to T?- . However 
if the number of C-side trees is odd the remaining dia- 
grams contribute to T//' . Hence it contributes to T,^. 

Taking into account the diagram connecting lead i and 
lead j directly, which has no Andreev reflection this dia- 
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gram contributes to Tfj^ and therefore to +2^^- Accord- 
ing to the diagrammatic rules this contribution is simply 
given by NiNj /[N-[<i{l + x)]. The transmission difference 
Tf- therefore read 



OO 

1=2 



(22) 



1=2 



with Af and Bf given by (|17ap and (|17b|) . respectively, 
and the same equation holds for T/j with an exchange 
of e and h yielding a replacement oi (f) by —((> in the 
coefficients A and B including an exchange P'^ O P^. 
If the numbers of channels of the superconductors are 
the same this also implies a symmetry in the exchange 
of electrons and holes and T?- — —T!) which reduces the 
number of variables by a factor of 2 and therefore makes 
the arising equations easier to solve. 

Since we consider the chemical potential of the super- 
conductors being the same as that of the second lead the 
contribution to the conductance is given by the transmis- 
sion coefficient for reflecting an electron entering the cav- 
ity from lead 1 back into lead 1 again. The conductance 
g, defined by Ii = eg {fj,i — /ig) /ttH, therefore reads 



27Vi / de l-f IP 



N{1- A + B) 



dl 
de' 



(23) 



Note that this setup induces an asymmetry due to the 
channel numbers: Since the chemical potential of the su- 
perconductor is the same as that of lead 2 one cannot 
exchange the two channels and therefore the solution will 
not (in general) be symmetric under the exchange of A^i 
and N2- 



A. Low temperatures 

1. One superconductor 

Let US consider first the simplest case of no phase dif- 
ference and the absence of magnetic fields. Moreover 
we consider sufficiently low temperatures to approximate 
—df/de « (5(e). The contributions of the side trees are 
therefore given by (fT^ . and we obtain 

Ni {2N2X^ + ANix^ + 8N2x'^ + k) 



2{i2x + 1) [Ni + N2) 
Ni {x + 1) K 



2y/l + 6x + x^ {N1+N2) 



(24) 



with K = N2x'^ + 2xNi + QxN2 + 2Ni + N2. 

We can easily compare our result for the conductance 
correction 6g = g—gd to the limiting case found for small 
superconductors in Ref. [13 by expanding our result in a 
Taylor series in a; = iVs/iVN. The first non- vanishing 
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FIG. 16. (a) The quantum correction to the conductance of 
an Andreev billiard with one superconducting and two nor- 
mal conducting leads (see Fig. 1(d) with the superconducting 
chemical potential being the same as that of normal lead 2) as 
a function of the number of channels of the superconductor, 
X = As /An, for N2/N1 -> (space dashed line), iV2/Ai = 0.2 
(dotted line), iV2/A^i = 1 (solid line), A2/A1 = 7.2 (dashed 
line) and N2/N1 — > 00 (dashed dotted line), (b) The size of 
the superconductor for which the conductance becomes min- 
imal as a function of the ratio N2/N1. 



term in the Taylor expansion is exactly the contribution 
found by Jacquod and Whitney^^: 



.(2) - 



Ni {N2 - 4,Ni) 



(25) 



Therefore for small superconducting leads the correction 
becomes negative if lead 1 carries a sufficiently large 
amount of modes compared to lead 2. The ratio N2/N1, 
for which the conductance has no longer a minimum and 
hence the conductance correction becomes no longer neg- 
ative for small x, is approximately N2/N1 w 4 as found 
by Whitney and Jacquod^S. Interestingly, Fig. 16( 
shows moreover a second change in sign of Sg, namely 
when increasing x, for fixed A^i > A^2/4. This charac- 
teristic crossover is beyond the treatment of Ref. 
arises from higher order diagrams as calculated here. The 
size of the superconductor for which the conductance 
becomes minimal is shown as a function of N2/N1 in 
Fig. 16(b) Moreover we again found a doubling of the 
conductance for Ng/Njq 00, i.e. g — 2Ni, indepen- 
dently of the ratio N1/N2 in alignment with previous re- 
sults for quantum dots with only one normal conducting 
leadSi. 



2. Two superconductors with a phase difference 

The effect of a phase difference between two super- 
conducting leads depends sensitively on the ratios x and 
N2/N1. The result is fairly similar to the phase depen- 
dence of the conductance of an normal conducting region 
with one normal conducting lead and two superconduct- 
ing leads with a phase difference (j) found by different 
approaches^*'^^. While for most combinations the effect 
of the superconductor decreases with increasing phase 
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FIG. 17. Phase dependence of the conductance correction, 
(a) X = 0.5 and N2/N1 -> (space dashed hne), iV2/iVi = 0.2 
(dotted hne), iV2/A^i = 1 (sohd hne), 7V2/A''i = 2 (dashed 
hne) and N2/N1 = 7 (dashed dotted hne). (b) iV2/iVi = 0.2 
and X = 0.1 (space dashed hne), x — 1 (dotted hne), x — 1.2 
(sohd hne), x = 1.5 (dashed hne) and x — 2 (dashed dotted 
hne). 



difference due to destructive interference, in some cases 
the conductance becomes a non-monotonic function be- 
tween = and = TT as shown in Fig. [T71 The phase 
difference may even cause a change of the sign of the 
conductance correction. In Fig. [17] this can be seen for 
the case iV2 = A^i and x — 0.5 as well as for the cases 
N2 = 0.2Ni and a; = 1 or a; = 1.2. 

However if the number of channels of the supercon- 
ducting leads are equal all possible combinations of x and 
N2/N1 have in common that the conductance correction 
becomes zero for a phase difference <j) — it and that the 
conductance is symmetric and periodic in (j) with period 
2tt. This can also be seen from p7a[ b) and (jH]): The sym- 
metry follows from the symmetry of Bi and P as well as 
the fact that Ai does not depend on explicitly. If (j) is 
replaced by -f- 27r the side tree contribution changes its 
sign. Ai and the first part of Bi as well as the last two 
terms of Bi for even Vs are symmetric in P and there- 
fore give the same contribution as before. However for 
an odd I the last two parts of Bi are antisymmetric in P 
but when increasing the phase difference by 27r the phase 
factors contribute an additional minus sign such that the 
total contribution stays the same. 

This symmetry towards (f> still holds if the two super- 
conducting leads provide different numbers of channels. 
Wc show the crossover from a setup with two supercon- 
ducting leads having the same number of channels to the 
setup with just one superconducting lead for different 
cases in Fig.[THl As in (|17bp y = {Ns^ - A^S2)/^S is the 
difference between the numbers of channels in the super- 
conducting leads such that y — corresponds to the sym- 
metric case A^Si = -^82 ^'^d y — ±1 to the case of just one 
superconducting lead. We found that in the cases where 
the conductance correction had a dip at cf) — ±7r first of 
all this dip vanishes and after that the conductance cor- 
rection tends to the (f) = result for all </) if y is increased. 
If there is no dip the result converges monotonically to 
the result for the case with just one superconducting lead. 
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FIG. 18. Dependence on the difference between the numbers 
of channels of the superconductors with x = 0.5 and y = 1 
(space dashed line), y = 0.95 (dotted line), y = 0.9 (solid 
hne), y = 0.5 (dashed line) and y = (dashed dotted hne). 
(a) N2/N1 = 0.2, (b) N2/N1 = 1, (d) N2/Ni = 2 (d) N2/N1 = 
7. 
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FIG. 19. Magnetic field dependence of the conductance cor- 
rection for (a) X = 0.1 and (b) x = 1 with N2/N1 — 0.2 (space 
dashed line), N2/N1 = 1 (dotted line, N2/N1 = 3 (solid hne), 
N2/N1 = 7 (dashed line) and N2/N1 -5> 00 (dashed dotted 
line) 



3. Weak magnetic field 

If a magnetic field is applied time reversal symmetry is 
broken. Therefore building side trees becomes less likely 
and their contribution vanishes as 6 cxd as can be 
seen in Fig. [191 In particular for small superconducting 
channel numbers Whitney and JacquodiS predicted that 
the conductance correction decreases at zero temperature 
with (1 -|- b^)~^. However, we find a nonmonotonic be- 
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haviour for the dependence of the conductance correction 
on the magnetic field similar to that found for the case of 
one normal conducting lead^ and for the magnetic field 
dependence of the excess current^. It may even happen 
that the conductance correction is negative for 6 7^ al- 
though it is positive at 6 = 0. Since the diagrammatic 
rules depend on rather than 6 the conductance is a 
symmetric function of h and therefore satisfies the On- 
sager relation^i^ for a two terminal setup g{h) = g(—b). 
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B. Temperature dependence 

To include the effect of finite temperature we have to 
include the energy dependence of the side tree contribu- 
tion and the central encounters. For zero phase the sixth 
order equation for the side tree contribution factorizes 
and one has to solve a quartic equation. At leading or- 
der in Ns/N^^ the temperature dependence is given by 
the generalised zeta function. However as can be seen 
from Fig. [20] when including higher order terms strong 
derivations from this behaviour may occur again. 

The integral in (|2T|) of course can not be performed 
analytically. Therefore we calculated the integral numer- 
ically using gaussian quadrature with a total accuracy 
of 10~^°. The integral has been truncated at e = 1000 
where 9 — k^T / E^: is again the temperature measured 
in units of the Thouless energy. In Fig. [^D] we plotted the 
conductance correction versus the temperature. 

As one might expect from the Sommerfeld expansion 
the conductance correction has a local extremum at 9 — Q 
but not necessarily a global one as can be seen from the 
sohd and the dash dotted curves in Fig. 20(a) For a cer- 
tain range of combinations of the ratios x and N2/N1 the 
conductance correction increases with increasing temper- 
ature although it is positive for T = 0. A similar effect 
known as the reentrance of the metallic conductance has 
previously found in NS-structures^i^. Moreover it may 
happen that the conductance correction changes its sign 
when the temperature is increased such as in the case 
N2/N1 — 0.2 and x = 1.4 shown by the solid line in 



Fig. 20(a) 



When including a phase difference again this affects the 
temperature dependence. It may again cause a nonmono- 
tonic temperature dependence even if it is monotonic at 
(j> = 0. Moreover in contrast to the case (p = where the 
sign of the conductance correction with increasing tem- 
perature may only change from positive to negative it 
can be the other way around for (f> ^ 0. The temperature 
dependence of the conductance correction for different 
phases is shown in Fig. [20r b).(c). 

In all cases for large temperatures the conductance cor- 
rection tends to zero. However this limit is approached 
only very slowly. 
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FIG. 20. The conductance correction as a function of the 
temperature Q = k^TjET- (a)iV2/iVi = 0.2, = and x = 
1.2 (space dashed line), x — 1.3 (dotted line), x = 1.4 (solid 
line), X = 1.5 (dashed line) and x = 1.6 (dashed dotted line), 
(b) (f) — n/i (space dashed line), (j) — 7n/20 (dotted line), 
(j) = ■k/2 (solid line), (j) = 137r/2Q (dashed line) and <j) = 3tt/4 
(dashed dotted line), N2/N1 = 1 and x = 0.5. (c) N2/N1 = 
0.2 and x = 1.3. 



VII. 



CONDUCTANCE OF A DOUBLE DOT 
SETUP 



For the double dot model shown in Fig. 1(b) [ where 
the quasiparticles stay on average a time St ■ td in the 
neck connecting the two dots, the necessary modifications 
are more substantial than in the previous case. First of 
all we have to find a way to calculate the transmission 
probabilities at all orders in x and in n = TVn/A^N where 
Nn is the number of channels of the neck connecting the 
two different dots. For simplicity we additionally assume 
that the two dots have the same dwell time. Since we 
consider the neck to be represented by an ideal lead every 
electron in the left dot entering the neck at channel i 
leaves the neck at the same channel into the right dot 
and vice versa such that a traversal of the neck yields a 
factor Nn- 

The idea of the calculation is as follows: We first of 
all split the whole setup into two distinct parts each hav- 
ing two normal leads (one of them being the neck) and 
one superconducting lead as it has been done in Ref. Issl . 
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FIG. 21. (a) An 3-encounter of the backbone touching the 
neck, (b) An 3-encounter of a side tree touching the neck. 



To find a recursion relation for the transmission we cut 
the backbone not only at the first encounter but also at 
the first traversal of the neck, i.e. the diagram is cut 
when the first path pair enters an encounter or traverses 
the neck without having entered an encounter. Since 
traversing the neck contributes a factor N^^ we will as- 
sign this factor to a path pair hitting the neck. Due to 
the splitting of diagrams at the neck we thus get path 
pairs starting at a neck which do not contribute a factor 
iVn. Thus a diagram starting at the neck in dot j con- 
tributes to T"-^/Nj. In order to simplify the handling 
of the contributed numbers of channels here we define a 
'normalised' transmission coefficient i"^^ — T°!j^ /Nj. 



Consider a diagram starting in lead 1 with an electron, 
thus contributing to with j £ {1, 2} and a G {e, h}. 
If the first path pair enters an encounter without having 
traversed the neck the diagrams are treated in the same 
way as above: The remaining part after cutting the dia- 
gram right after the first encounter contributes to if 
the number of C-side trees emerging from the encounter 
is even and to t"^ otherwise. If the number of ^-side trees 
is even the contribution of the first path pair, the first 



encounter and its side trees is given by 



(26) 



where the subscript '11' indicates that the encounter lies 
fully in dot 1. P" , i £ {1, 2}, a G {e, h}, is the contribu- 
tion of side trees starting with an a-type quasi particle 
in dot i and will be derived below. 

If the e-e* path pair leaving the first encounter and the 
odd numbered side trees attached to the first encounter 
traverse the neck before entering another encounter the 
first encounter can touch the neck. As shown in Fig. 21(a) 
the odd numbered side trees and the e-e* path pair leav- 
ing the encounter then start in dot 2, since we treat the 
neck as a ballistic connection between the two dots, while 
the total number of side trees stays the same. Hence the 
contribution becomes 



^21,1 — 



l-l 



[l + i{2p-l + l)e] 



[P^e^' 



P^Y {P^e'^'^Pl^y 



(27) 



Since the e-e* path pair now starts in dot 2 the remaining 
part of the diagram in this case contributes to tf2- The 
additional phase factor e"^'''^ is the phase accumulated by 
an e-e* path pair which traverses the neck. 

For the encounters with an odd number of side trees 
built by C we have to be even more careful. This is 
because the encounter now can be moved into either 
the neck or into the superconductor coupled to the dot 
the encounter lies in. The fact that when an encounter 
touches the neck the odd numbered side trees have to 
traverse the neck means that the path pairs of the back- 
bone entering and leaving the first encounter have both 
to be in the same dot after traversing the encounter. The 
I — 2 in total side trees still starting in the left dot are the 
originally even numbered side trees and therefore start 
with a hole. 

If again the first encounter is in the left dot the con- 
tribution of the encounters with an odd number of side 
trees built by ( is 



Bki = - E [1 + i (2P - ^ + 2) e] {P^r' {Plf [{Pi) 

P=0 [ 



l-p-2 



Si 



l-p-2 



^.[2p-l+2)e5r ^phy ^phy (P|)^+^ [(P2')*] 



l-p-2 



1 l~P~l 



(28) 
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FIG. 22. A side tree starting in dot 1 may also have parts lying in dot 2. If a quasiparticle hits the neck we will treat this 
as an electron being retroreflected as its counter-particle indicated by the empty boxes and include an additional factor of the 
side tree starting in the other dot for each traversal of the neck, 'e' and 'h' denote that the side tree starts with an electron or 
hole, respectively. 



Similar formulae as above hold for the remaining coefficients Afj j and B"j j, i,j € {1, 2}, a G {e, h}. 

If the e-e* path pair hovifever hits the neck without having traversed an encounter we cut the diagram at the 
neck such that the first part contributes Nn/ {Ni + Ng^ + N^) and the second part contributes to Therefore the 
'normalised' transmission coefficient is given by 




B,.it?r (29) 



1=0 



where j = 3 — j as well as /3 = h ii /3 — e and vice versa. 

This gives in general a 16 dimensional system of lin- 
ear equations which decomposes into four independent 
systems of linear equations. 

The complication of treating the double dot setup is 
that the side tree contributions are also different. We 
may calculate their contribution though following the 
steps in Sect. IIIII with slight changes. First of all it is 
no longer enough just to consider the generating func- 
tions / and / for side trees starting with an electron or a 
hole, respectively, but we have to consider the generating 
functions /i , /i , /2 and /2 for side trees starting with an 
electron (without the hat) or a hole (with the hat) in the 
left (subscript '1') or in the right dot (subscript '2'). Here 
we will only consider /i for the side trees starting with an 
electron in the left dot explicitly since the derivation of 
the remaining generating functions is similar and needs 
only simple replacements. 

First of all we have to slightly modify the way we look 
at the side trees. In general they may consist of parts 
lying in dot 1 and parts lying in dot 2 as indicated in 
Fig. [211 By cutting a tree starting in dot 1 at the neck 
we find that an electron traversing the neck builds a tree 
starting in dot 2 and finally comes back as a hole since the 
number of Andreev reflections provided by a tree is odd. 
Thus every time a link of a tree starting in dot 1 hits the 
neck we cut the tree at the corresponding channel and 
insert a retro-refiection there contributing an additional 
factor given by the side tree starting in dot 2, hence a 
factor /2 if an electron hits the neck and /2 if a hole hits 
the neck. Moreover we have to include the additional 
phase eSr due to the time a quasiparticle spends in the 



neck. Therefore we may write a diagrammatic rule for a 
path pair in dot 1 hitting the neck: 

• An electron hitting the neck contributes a factor 
7V„e'^^V2. 



• A hole hitting the neck contributes a factor 

f2. 



We may then start with (j6ap and adapt it to our new 
problem since the steps leading to this can be used here in 
exactly the same way. The second term in the first line 
and the last term in the second line which previously 
corresponded to a tree whose first encounter hits S2 is 
now given by a tree with its first encounter hitting the 
neck and therefore building a side tree in dot 2 yielding 
a contribution /2 . Therefore we replace the contribution 



(2) 



of an ^-encounter touching S2 by the contribution 
of an ^-encounter touching the neck. Thus 



fi = - ii^e-'-^/^ + J^e''^^ U 



1=2 



(30) 



where 7V(*) = Ni + TVg. + Nn, i e {1, 2}. 

Due to the assumption that the neck may be repre- 
sented by an ideal lead an Z-encounter may touch the 
neck even if the odd numbered side trees do not have 
zero characteristic. The only restriction to them is that 
they traverse the neck before having an encounter or hit- 
ting a superconductor. Thus when sliding an /-encounter 
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into the neck we get the situation depicted in Fig. 21(b) 
The odd numbered side trees, which start with an elec- 
tron, now start in dot 2 instead of dot 1. Moreover there 
are / path pairs traversing the neck each giving a phase 
e5T. Therefore we get 



Jn,l) _ 



or /• 

/: 



(31) 



xij and z'^J'j are obtained in the same way as in Sect. IIIII 
and are given by 



^o,l 



7V( 



(1 + i;e)f' 
^(_i)'e""*/2f'-i. 



(32a) 
(32b) 



The remaining steps then are again the same as in 
Sect, mil For a side tree starting with an electron in dot 
2 we have to exchange the labels '1' and '2' as well as 
the phase by —cj). If we consider side trees starting 
with a hole in dot i instead of electrons we have also to 
reverse the phase with respect to fi and exchange fj O 
fj, j S {1,2}. All in all after performing the sums using 
geometric series, making the change of variables ([5]) and 
setting r = 1 to get the side tree contributions out of 
gi and gi, we obtain 

(^_^peph _ (P^pl'f + 2iePf - l) Pf Pj^ 



pephy 



7V(1) 



,i0/2 



— \e5r 



(33) 



and similar equations for P/\ P| and Pj*. 

For simplicity in the following we consider zero tem- 
perature, i.e. e = 0, and equal leads Ni — N2 — 
and TVsi = Ns^ = Ns/2. Then one finds that P^ ^ Pf 
where a labels a hole if a labels an electron and vice 
versa. Moreover in this case exchanging electrons and 
holes is the same as exchanging dot 1 and dot 2 and is 
related to the fact that when exchanging electrons and 
holes this is essentially an exchange oi (p ■n- —cj). This 
also reduces the number of linear equations for the nor- 
malised transmission coefhcients by a factor of 2 since it 
yields A? = and Bf, = Pi where i = 2 if i = 1 and 

vice versa. Therefore we have III = B^^ etc. The set of 
4 nonlinear equations indicated by (j33|) decomposes into 
two copies of sets of two nonlinear equations where the 
second set is the same as the first one. 

The classical conductance can be found by using a sim- 
ilar recursion relation for the 'normalised' transmission 
coefficient as above which is for the diagonal diagrams 
given by 
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FIG. 23. The conductance correction of the dobule dot setup 
wit (a) n = 0.5 and x — 5 (dashed line), x — 20 (solid line) 
and X = 100 (dotted line), (b) n — 20 and x = 1 (dashed line), 
X — 20 (solid line) and x = 100 (dotted line), (c) x = 0.2 and 
n — 0.3 (dotted line), n = 0.5 (solid line) and n = 0.7 (dashed 
line), (d) n — 5 (dotted line), n = 20 (solid line) and n = 100 
(dashed line). 



With the classical conductance obtained using this recur- 
sion, gel = N^{2n + x)/{2 + Ax + 8n), where n = iV„/iVN, 
X = Ns/N-^ and iVs = Ns^ + Ns2, at zero temperature 
we get the conductance correction shown in Fig. 1231 The 
conductance correction is symmetric in the phase differ- 
ence due to the symmetry in exchanging the two dots and 
the symmetry with respect to the exchange of electrons 
and holes. Moreover we found that the conductance cor- 
rection is periodic in the phase difference with a period 
of 27r as it can be seen from and ([55]) : If we re- 

place <^ — > 4> + 2Tr the side tree contributions change their 
sign. However in Af^ ; and the first and the last term 
of Bfj I in (pS)) this sign cancels since the total number 
of side tree contributions entering them is always even. 
An additional minus sign enters the second term of P° ^ 
if I is odd, but when I is odd the phase factor of this 
term also contributes an additional minus sign. If I is 
even the total number of side trees entering this term is 
again even and the phase factor does not contribute an 
additional minus sign on increasing the phase difference 
by 2n. Thus the transmission coefficients are symmetric 
under replacing (f) hy (f) + 2Tr . 

Moreover we find that for a small number of channels 
in the neck the conductance correction is negative. It 
(34) increases if n is increased and finally becomes positive. 
For such values of n where the conductance correction 



19 



is positive we also find that the conductance correction 
vanishes for (j) = ±7r. For smaller values of n the changes 
in the conductance correction when changing a; are most 
pronounced around ±7r as it can be seen from figures 
[2li,b. 



VIII. THERMOPOWER OF THE SYMMETRIC 
AND ASYMMETRIC HOUSE 



the numbers of channels of the two normal leads be- 
ing equal, since these numbers only enter by a prefactor 
N1N2/ {Ni + N2). It is further assumed that no magnetic 
field is applied. Note that we consider the superconduc- 
tors to be isolated thus adjusting their chemical potential 
such that the net current through the superconducting 
leads is zero. 



If the normal leads additionally have different temper- 
atures there is also a coupling between the electrical cur- 
rent and the temperature difference. For a two terminal 
setup with the superconductors being isolated and the 
two normal leads having a temperature difference AT 
and a voltage difference AV the electrical and thermal 
current are given by2S 



G B 

r s 



Ay 

AT 



(35) 



where G = 2e^g/h is again the conductance essentially 
given by p9| and 



(36a) 




(36b) 



(36c) 



where T 



- r de{df/de)Tf and T" 

- de e {df/de) T^f . Note that in Ref.lH these coeffi- 
cients have been written down only for low temperatures 
using the Sommerfeld expansion. 

An estimate of the thermo-electric coupling is for ex- 
ample provided by the thermopower 



1=0 



AV 
AT 



7=0 



B 
G' 



(37) 



The calculation of the thermopower is therefore closely 
related to the electrical transport since due to ([57]) we 
have to evaluate the same transmission coefficients. We 



will first consider the two cases shown in Fig. 1(e) and 



which are topological equivalent to the symmet- 
1(e) I and asymmetric house (Fig. 1(f) )^2^. We 



Fig.|l(f) 
ric (Fig. 

will restrict ourselves to Andreev billiards with two nor- 
mal leads and two isolated superconducting leads with 
a phase difference (j). The results presented here are 
again in leading order in the channel numbers and for 



A. Symmetric house 



We start with the setup Jacquod and Whitney called 
the symmetric house^^ shown in Fig. 1(e) They treated 



the transmission coefficients perturbatively in the ratio 
X = Ns/N]S! up to second order. Within this approxima- 
tion they argued that for the symmetric house with equal 
numbers of superconducting channels the thermopower 
vanishes in second order in x since B, defined in (I36a|) . is 
antisymmetric in exchanging electrons and holes which 
yields an exchange e ^ — e and in leading order in 
N is equivalent to reversing the superconducting phase 
(j> — ^ —(p. Since additionally the electrical conductance G 
is symmetric as already seen in Sect. |V] the thermopower 
S — —B/G is antisymmetric in the phase difference. On 
the other hand the result has to be symmetric under ex- 
changing the superconducting leads which is again equiv- 
alent to reversing the sign of (j). Thus the thermopower 
has to be zero. With our approach we find that this ar- 
gument holds to all orders and may also be seen from the 
discussion in Sect. IIVI There we stated that as long as 
the superconducting leads both provide the same number 
of channels the transmission probability is symmetric in 
exchanging electrons and holes which when inserted in 
(PSa|l gives B = 0. 

However if iVg^ ^ Ng^ the symmetry under exchanging 
the two superconductors is broken and the thermopower 
does not vanish. By using the transmission coefficients 
found in Sect. IIVI inserting into (IT51) as already done 
in Sect. |V] and (j36ap and performing the integrals nu- 
merically with a total accuracy of 10~^'^ we find the re- 
sults shown in Fig. 1241 The thermopower is found to 
be antisymmetric and 27r-periodic in the phase differ- 
ence (j). When the difference y = (A^Si — -^S2)/^s is 
increased starting at zero a non-zero thermopower ap- 
pears which increases up to about y — 0.5. If the differ- 
ence is increased further the resulting thermopower de- 
creases again and vanishes at y = 1 which corresponds to 
the case of just one superconductor. Increasing the total 
number of superconducting channels x = Ng /A^n results 
in an increase of the thermopower while an increase in 
the temperature 6 = kBT/ET causes a decrease of S. 
This is due to the fact that the thermopower of Andreev 
billiards is a pure quantum mechanical property since the 
diagonal approximation gives S = 0. 
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FIG. 24. The thermopower of the symmetric house in units 
of 10~"^fcB/e with different numbers of superconducting chan- 
nels: (a) Dependence of the thermopower on the difference be- 
tween the numbers of superconducting channels with x = 0.2, 
6 = 0.2 and y = 0.3 (dashed line), y = 0.5 (solid line) and 
y = 0.9 (dotted line), (b) Dependence on x with y = 0.3, 
6 — 0.2 and x = 0.2 (dotted line), x = 0.5 (solid line) and 
X — 0.7 (dashed line), (c) Dependence on the temperature 
with X = 0.2, y = 0.3 and 6 = 0.2 (dashed hue), 6 = 1.5 
(solid line) and — 2 (dotted line). 



B. Asymmetric house 

Moreover for N^^ = iVs^ we can also generate a non- 
zero thermopower by inserting a neck at one of the two 
superconducting leads, say at Si as in Fig. 1(f) in which 
the trajectories spend an additional time St ■ t^. An 
e-h pair hitting the superconductor Si picks up an addi- 
tional phase cSt. Thus the total phase provided by the 
neck plus the Andreev reflection at Si is (— 0-1- 2e(5r)/2 if 
an electron is converted into a hole and {(j> -\- 2£St)/2 if a 
hole is converted into an electron. Therefore the electron 
hole symmetry is broken leading to an asymmetry in the 

phase difference as well as in the energy. So T^'^ ^ T^j'^ 
with a (/3) labelling a hole if a (/3) labels an electron and 
vice versa. Note that we treat the neck as an ideal lead 



such that every quasiparticle entering the neck hits the 
superconductor before leaving the neck again^. When 
redoing the steps of sections IIIII and IIVI for the calcu- 
lation of the transmission coefficients again one finally 
finds the following changes: The variable /3 = cos{(f>/2) 
has to replaced by P'^ — cos[(0 — eST)/2] for a side tree 
starting with an electron and by = cos[(0 -I- e(5r)/2] 
for a side tree starting with a hole. Furthermore the side 
trees starting with an electron are multiplied by e~"^*'^/^ 
and those starting with a hole by e"^"^"^/^. However since 
in Af and B" each factor P" is paired either with a fac- 
tor or with a factor (P")* these additional factors 
cancel. Additionally if the incident quasiparticle is an 
electron the phase (f> in the second term of (jl4bp also has 
to replaced hy 4> — 2eST. Again for an incident hole the 
phase (p has to be replaced hy —cj). 

All in all we have to solve a 4 x 4 system of linear 
equations rather than a 2 x 2 system of linear equations 
as it was for the conductance of the symmetric version. 
Using Ni = N2 = AfN/2 the four equations are 
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with Ai being the same as in (jl4ap but with P'^ and 
P'* redefined to depend on Z?*^ and /?'', respectively, and 
6 = 0. Analogously one finds 
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I 

By inserting the transmission coefficients into (fT9|) and ically using Gaussian quadrature with a total accuracy 
P6ap we find the thermopower S — —B/G. To this end of 10~^°. The results for the thermopower for different 
we again integrated the transmission coefficients numer- values of 5t, x and different temperatures is shown in 
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FIG. 25. The thermopower in units of 10 ^ks/e as a function 
of the phase difference 4>. (a) St = 0.03, 6 = KbT/Et = 0.2 
and X = 0.2 (dotted hne), x = 2 (solid Une) and x = 10 
(dashed hne). (b) St = 0.03, x = 0.2 and 6 = 0.2 (dotted 
hne), 6 = 4 (sohd hne) and 9 = 10 (dashed hne). (c) 9 = 0.2, 
X = 0.2 and St = 0.03 (dotted hne), St = 0.1 (sohd hne) and 
St = 0.3 (dashed hne). 



Fig.m 

We find that the antisymmetry in phase found by 
Whitney and Jacquod in second order in x holds up to 
all orders in x. However this antisymmetry is in contra- 
diction to previous experimental measurements^ii^ on 
diffusive normal metal regions coupled to two supercon- 
ductors. Moreover one can see from Fig. [55] that the 
thermopower is 27r-periodic in (jj. The arguments for the 
periodicity of the conductance in Sect. |V] also apply to 
B and thus the thermopower is periodic in the phase 
difference <j) with period 27r. A period of 27r had also 
been found previously in Ref. Eol . The antisymmetry and 
the periodicity may also be obtained by (I39a[ b) com- 
bined with p9)) and (|36a[) : Due to the summation over 
p, Bf and B^ are symmetric under an simultaneous ex- 



change (j) O —0 and e ■H' —e and satisfy Bj^ = Bf 
Thus replacing by —0 is the same as replacing e by 
— e. Now if we replace e by — e we get an additional mi- 
nus sign in p6ap and therefore a minus sign in the ther- 
mopower. In contrast to the symmetric case however the 
symmetry due to the exchange of the leads is broken and 
thus the thermopower is non-zero but antisymmetric in 
(j). Furthermore for specific combinations of x and St the 
thermopower as a function of the phase difference may 
show additional oscillations with period smaller than 2tt 
(see Fig. 25(a) |. However these additional oscillations 



are smoothed out if the temperature is increased. 

Over a pretty wide range of increasing the temper- 
ature the thermopower also increases before it shrinks 
again. This increase of the thermopower with increasing 
temperature is related to the decrease of the thermal re- 
sistance found in Ref. [s^ Fig. 25(c) also shows that the 



thermopower increases if the time a quasiparticle spends 
in the neck is increased. This is due to the fact that with 
increasing St the electron-hole symmetry is increasingly 
broken. 

It should be noted that the thermopower again arises 
solely from the non-diagonal diagrams since for the di- 
agonal diagrams the energy differences of the links are 
always zero and thus the neck does not play any role. 



IX. CONCLUSIONS 

We have shown in this article that the transport prop- 
erties of chaotic cavities are strongly affected by the in- 
troduction of superconducting leads. In particular, using 
a semiclassical framework, we have extended the work of 
Ref. [13 to all orders in the ratio Ns/Nys of the sizes of the 
superconducting and normal leads. With this we could 
reproduce the large- iV limit of the random matrix the- 
ory result^^ for two isolated superconductors which shows 
in the limit Ng/N-^ oo the conductance doubling of 
the N-S interface consisting of one normal conducting 
and one superconducting leadi. Our result shows a 27r- 
periodicity and symmetry with respect to the phase dif- 
ference of the superconductors which have already been 
found by other approache a'^'^^ . If a magnetic field is ap- 
plied or the temperature is increased the conductance 
correction will decrease towards zero. 

An Andreev billiard with two normal leads and two 
superconducting leads having the same chemical poten- 
tial as one of the two normal leads shows more interest- 
ing features than the setup with superconducting islands. 
Depending on the numbers of channels, these may be sim- 
ilar to those already found previously in the conductance 
of N-S junctions. Again the conductance is doubled in 
the limit Ng /Njq — >■ oo and is 27r-periodic and symmetric 
with respect to the phase difference which has been ob- 
served for N-S junctions with two superconducting leads 
in several approachea^SiS^ii^ii^. However the conductance 
shows a non-monotonic behavior similar to the conduc- 
tance through quantum dots with one normal and two 
superconducting leads^S,. The phase difference may even 
cause a change of the sign of the conductance correction. 
Furthermore the magnetic field dependence also inher- 
its a non-monotonic behavior from the N-S junction^ii^ 
and again the sign of the conductance correction may 
be changed by increasing the magnetic field. This has 
already been found by Whitney and Jacquod^'' in their 
consideration of the contribution up to second order in 
A's/A'n. a non zero temperature may even cause an in- 
crease of the conductance correction similar to the reen- 
trance in the case of the N-S junction^i^. 

Separating the Andreev billiard into two dots each con- 
sisting of one normal and one superconducting lead and 
connecting them by a neck necessitates a splitting of the 
diagrams when traversing the neck. A small number of 
channels in the neck causes a reduction of the conduc- 
tance while for a large neck the superconductors enhances 
the conductance as long as the phase difference between 
the two superconducting leads cf) ^ n. The conductance 
of this setup is again symmetric and periodic in (j) with a 
period of 27r. 

The investigations of the thermopower in Ref. have 
also been extended to all orders in the number of su- 
perconducting channels for the symmetric house and the 
asymmetric house. For these two setups we could show 
that the antisymmetry of the thermopower towards the 
phase difference (p holds in all orders in Ng/Nj^i and that 
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the thermopower of the symmetric house with equal num- 
bers of superconducting channels is identically 0. In the 
case of the asymmetric house this antisymmetry how- 
ever is in contradiction with some previous experimental 
results^iS^, which found a symmetric thermopower for 
diffusive normal regions. But the thermopower oscillates 
in both cases with a period 2tt in agreement with ex- 
perimental results^i^. Additionally the thermopower 
increases with increasing temperature over a pretty wide 
range. With the treatment represented in Sect. I Vlll where 
we calculated the conductance of the double dot we would 
in principle be able to calculate also the thermopower of 
this setup. However including the necessary energy dif- 
ference and the difference between the numbers of chan- 
nels of the two superconductors increases the complexity 
of the equations to such an extent that we have not got 
a reasonable solution of the four coupled nonlinear equa- 
tions for the side tree contributions. 

Like for the the density of states^^, the considerations 
in this article show that in the case of normal metal- 
superconductor hybrid systems, a semiclassical treat- 
ment based on the diagonal approximation is not suf- 
ficient to describe such systems. Instead one has to con- 
sider all diagrams consisting of path pairs and encoun- 
ters contributing in leading order in the channel number. 
Importantly it shows that the effect of superconductors 
on attached normal regions is of the order of the total 
number of channels and the leading order corrections in 
-^S /-^N compete with higher order terms for larger num- 
bers of superconducting channels. 

However up to now there are neither confirming nor 
disproving experimental data. Thus it would be in- 
teresting to see results of measurements upon the con- 
ductance of Andreev billiards. While the realization of 
ballistic Andreev systems based on InAs has been ex- 
perimentally shown^, corresponding transport measure- 



ments have yet to be performed. 

The results presented here are only valid in leading 
order in the inverse channel numbers, and further in- 
vestigations could include the calculation of subleading 
order terms, especially the weak localisation correction. 
Doing this first of all requires a way of systematically 
finding all the diagrams contributing in the next order in 
1/iV. The key problem then is to find a suitable recursion 
since the structures found here break down. Namely the 
trees have no longer to be trees and the diagrams may 
not necessarily have a backbone. However the types of 
recursion relations presented here would prove useful in 
such a case. We may also wonder about the effects of 
superconductors on the noise and other finer transport 
statistics. 

Here we considered Andreev billiards coupled to at 
most two superconductors. However extending our cal- 
culation to a higher number of superconductors seems 
(maybe up to solving the equation found for the side tree 
contributions) straight forward. We also restricted our- 
selves to zero Ehrenfest time. However it has been shown 
that the Ehrenfest time plays a crucial role in Andreev 
billiards^^. Allowing non-zero Ehrenfest time yields a 
simple replacement for the side treesi^'i^ but the effect 
on the backbone has to be investigated further. 
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